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Abstract 

In this article, we take the point of view that the light scalar meson ao(980) 
is a conventional qq state, and calculate the coupling constants g aoV7T o and 
9a r]'ir with the light-cone QCD sum rules. The central value of the coupling 
constant g aorj7T o is consistent with the one extracted from the radiative decay 
0(1020) — > ao(980)7 — > rjir ^. The central value and lower bound of the decay 
width T ao _ trjn o = 127^4g MeV are compatible with the experimental data of 
the total decay width r ao ( 980 ) = (50 — 100) MeV from the Particle Data Group 
with very model dependent estimation (the decay width can be much larger) , 
while the upper bound is too large. We give possible explanation for the 
discrepancy between the theoretical calculation and experimental data. 

PACS numbers: 12.38.Lg; 13.25.Jx; 14.40.Cs 
Key Words: a (980); Light-cone QCD sum rules 

1 Introduction 

The light flavor scalar mesons present a remarkable exception for the constituent 
quark models, the structures of those mesons have not been unambiguously deter- 
mined yet [U [2J [3j H] . Experimentally, the strong overlaps with each other and the 
broad widths (for the / (980), a (980), / (1710), the widths are relatively narrow) 
make their spectra cannot be approximated by the Breit-Wigner formula. The nu- 
merous candidates with the same quantum numbers J PC = ++ below 2 GeV can not 
be accommodated in one qq nonet, some are supposed to be glueballs, molecules and 
multiquark states [21 SI - The more elusive things are the constituent structures 
of the mesons /o(980) and ao(980) with almost the degenerate masses. 

In the naive quark model, a = (uu—dd) / v2 and fo = ss; while in the framework 
of the tetraquark models, the mesons /o(980) and a (980) could either be compact 
objects (i.e. nucleon-like bound states of quarks with the symbolic quark structures 
fo = ss(uu + dd) I y/2 and a = ss(uu — dd) / y/2 [SJ [6]) or spatially extended objects 
(i.e. deuteron-like bound states of hadrons: KK molecules [3 [8]). The hadronic 
dressing mechanism takes the point of view that the mesons /o(980) and ao(980) have 
small qq cores of typical qq meson size, strong couplings to the intermediate hadronic 
states (KK) enrich the pure qq states with other components and spend part (or 
most part) of their lifetime as virtual KK states [91 [101 [11]. In the hybrid model, 
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those mesons are tetraquark states {qq)i{qq)z in the S- wave near the center, with 
some constituents qq in the P-wave, but further out they rearrange into (<?<?) 1 (<?<?) 1 
states and finally as meson-meson states [21 Hj- All those interpretations have both 
outstanding advantages and obvious shortcomings in one or other ways. 

We can study the scalar mesons through their couplings to two pseudoscalar 
mesons, two-photon decays and radiative decays. The radiative decays 0(1020) — > 
7r°7r°7 and 0(1020) — > rjir ^ have been the subject of intense investigation [121 [EH 
dU [T51 [TBI El EE]- From the invariant 7r°7r° and f]ir mass distributions, we can 
obtain many information about the nature of the /o(980) and ao(980) respectively. 

In this article, we take the scalar mesons ao(980) and /o(980) as the conventional 
qq states, and calculate the values of the coupling constants g aorpr o and <7 O0 ^7r with 
the light-cone QCD sum rules. The coupling constant g aoV7T o is a basic parameter in 
studying the radiative decay 0(1020) — ► a (980)7 — ► r]7i°'j. In previous works, the 
mesons /o(980), ao(980), D s0 , D s i, B s0 and B s i were taken as the conventional qq, 
cs and bs states respectively, and the values of the coupling constants gf KK, ga KK, 
gD s0 DK, gD sl D*K, gB s0 BK and gB sl B*K have been calculated with the light-cone QCD 
sum rules [EH [201 [221 [231 [21] • The large values of the coupling constants support 
the hadronic dressing mechanism. In Ref . [25] . the authors study the coupling con- 
stant g aori TT with the interpolating current J® = 4= [u-y^j^u + d^^y^d — 2^7^755] , 
a complex subtraction procedure is taken due to the asymmetric Borel parameters 
Mf 7^ Mf . In this article, we study the coupling constants g aon!r o and g ao n'n° to- 
gether, and a simple subtraction procedure is taken. The decay /o(980) — > nn can't 
occur at the tree level if the scalar meson / (980) is a pure ss state, it should have 
some nn components, the coupling constant <7/ o7nr has also been calculated with the 
light-cone QCD sum rules [26J. 

The light-cone QCD sum rules approach carries out the operator product ex- 
pansion near the light-cone x 2 ~ instead of the short distance x ~ while the 
nonperturbative matrix elements are parameterized by the light-cone distribution 
amplitudes instead of the vacuum condensates [2Z1 [2S1 1221 EDI EI] • The nonpertur- 
bative parameters in the light-cone distribution amplitudes are calculated by the 
conventional QCD sum rules and the values are universal [321 E31 El]- 

The article is arranged as follows: in section 2, we obtain the coupling constants 
ga Q r]TT and g ao ri'n with the light-cone QCD sum rules; in section 3, numerical results; 
section 4 is reserved for conclusion. 
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2 Coupling constants g aQT]7l o and g aQr] ' n o with light- 
cone QCD sum rules 

In the following, we write down the definitions for the coupling constants <? ao777r o and 

9aor)'ir° j 



(a \r]7T ) = ig aoVn a = i\j^g , 
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(a \r]'7i°) = ig aWw o = i—g , (1) 

where we have used the phenomenological lagrangian C = gTr[SPP], the S and 
P stand for the light nonet scalar mesons and pseudoscalar mesons respectively. 
We study the coupling constants g aonK o and g am '^ with the two-point correlation 
function H^ip, q), 

U,(p,q) = i j ' d A xe-^^\T{J^)J{x)}\^{p)), (2) 
J^x) = u(x)^^u{x) + d(x)^^d{x) , 

j( x j = u(x)u(x) - d(x)d(x) ^ 

where the currents J^{x) and J(x) interpolate the pseudoscalar mesons rj, rj' and 
scalar meson ao(980), respectively, the external 7r° meson has four momentum p^ 
with p 2 = m 2 . One may think that it is more convenient to take the octet current 
J^(x) and singlet current J°(x) 

8 u(x)-f^ b u(x) + d(x)>y^d(x) - 2s(x)-/^ 5 s(x) 

jq = u(x)^ 5 u(x) + d(x)^^ b d(x) + sjx^^sjx) 

to interpolate the pseudoscalar mesons rj and rj' respectively. The ss components of 
the interpolating currents have no contributions at the level of quark-gluon degree's 
of freedom, the octet current J^(x) and singlet current J^(x) lead to the same 
analytical expressions. The J^{x) is a linear composition of the octet current J^(x) 
and singlet current J®(x), we choose it to interpolate the mesons rj and rj' together, 

= + . ( 5 ) 

Despite which interpolating currents one may choose, the couplings with the a (980)7r° 
take place through the uu and dd components of the pseudoscalar mesons i] and rj' 
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(not the ss component) at the level of quark-gluon degree's of freedom. Although 
the coupling constant g aov > w o has no direct phenomeno logical interest, we take into 
account the rf meson to facilitate subtractions of the continuum states and obtain 
more reliable QCD sum rules, we will revisit this subject at the end of this section. 
The correlation function ILj(p, g) can be decomposed as 



n M (p,g) = iU(p,q)q^ + iU A (p,q)p^ 



(6) 



due to Lorentz covariance, we choose the tensor structure for analysis. 

According to the basic assumption of quark-hadron duality in the QCD sum rules 
[521 [321 El], we can insert a complete sets of intermediate hadronic states with the 
same quantum numbers as the current operators J^{x) and J(x) into the correlation 
function II^p, g) to obtain the hadronic representation. After isolating the ground 
state contributions from the pole terms of the mesons rj, rf and ao(980), we get the 
following result (we present some technical details in the appendix), 



n, 



(v\J(o)\A + 



(0|^(0)|r7(?+P))/„,^M,ox , (0\M0)\y'(q + p)) 

M 2 - (q + p) 2 -ie 

i 



M% - {q + pf -ie 

i2f rj (q + p) ll 
V6[M 2 -(q + p) 2 -ie\ 

i2f rj ,(q + p) ll 
V^[M 2 -(q + p) 2 -te] 

i2gf v f ao M ao 



<t/|J(0)|7r°> 



(r]\a TT ) 



q 2 -M 2 +te 



(r]'\a ir )- 



(ao(?)|J(0)|0> + 
(a (g)|J(0)|0) + 



3[Af2-( 9 + p) 2 -ie] 



at) 



(P + ?)/* + 



3[M*-(q+ P y-ie] [M^-qt-ie] 

where the following definitions for the weak decay constants have been used, 

«2/„ 



(7) 



(o\Mo)\v(p)) 
(W)W(p)) 

(0\J(0)\a o (p)) 



i2f rj 



Pn, 
Pn 



fa M aQ . 



(8) 



I uu+dd— 2ss 



), rf 



We have take the ideal mixing limit for the rj and r]'(i.e. rj 

\ uu+ ^ +ss )), and neglect the anomaly contribution. 

In the following, we briefly outline the operator product expansion for the corre- 
lation function iT M (p, q) in perturbative QCD theory. The calculations are performed 
at the large space-like momentum regions (q + p) 2 <C and g 2 C 0, which corre- 
spond to the small light-cone distance x 2 ~ required by the validity of the operator 
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product expansion approach. We write down the propagator of a massive quark in 
the external gluon field in the Fock-Schwinger gauge firstly [35J, 



1 



' U ;5ij - I dvg s G^iyxx + (1 - v)x 2 ) 



k 2 -m 2 lJ / ""^'"w 



o 



"1 M + m 1 \ 1 ] /^\ 

Substituting the u and d quark propagators and the corresponding 7r-meson light- 
cone distribution amplitudes into the correlation function Ti^p, q), and completing 
the integrals over the variables x and k, finally we obtain the analytical expressions. 
In calculation, the two-particle and three-particle 7r-meson light-cone distribution 
amplitudes have been used [361 S3 SHI EH], the explicit expressions are given in 
the appendix. The parameters in the light-cone distribution amplitudes are scale 
dependent and are estimated with the QCD sum rules [361 Ell EE EH] . I n this article, 
the energy scale /i is chosen to be \i = 1 GeV. 

After straightforward calculations, we obtain the final expression of the double 
Borel transformed correlation function II at the level of quark-gluon degrees of free- 
dom. The masses of the pseudoscalar meson and scalar meson are M v i = 0.958 GeV 
and M ao = 0.985 GeV respectively, 

M 2 , 

~ 0.49 , (10) 



M-, + M 2 

r)' 1 ao 



there exists an overlapping working window for the two Borel parameters Mf and 

Ml Ml 
mJ+m% ■ 



Mf , it is convenient to take the value Mf = Mf , M 2 = -^tt^i- We introduce the 



threshold parameter sq and make the simple replacement, 



to subtract the contributions from the high resonances and continuum states [35J 
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Finally we obtain the sum rule for the coupling constant g, 




M 2 

$(1 — a — a g , a g , a) 



where 



Uq 



A\\(ati) + A±(ai) - V\\(oii) 
m « + M o(l -u )ml, 

Ml 
Ml + Mf ' 



V\ (a; 



(12) 



and we have taken the isospin limit m u = ra^. 

In Ref. [25] (also in Refs. [22l [23l [26]). a complex subtraction procedure is taken 
due to the asymmetry Borel parameters, Ml ^ Mf. In the light-cone QCD sum 
rules, we often take the technique developed in Ref. [35] to obtain the spectral den- 
sities at the level of quark-gluon degrees of freedom, 



n 



o A - (q + up) 2 
Pqcd(s 



du 



7 ds 



A [s ~ (p + q) 2 ] [s - q 2 ] 

f°° PQCpfcl, 32)5(3! - Sg) 

Ai Ja 2 [*i - (p + g) 2 ] [s 2 - <? 2 



■ds\ds2 



(13) 



where the /(tt) are functions of the two-particle light-cone distribution amplitudes, 
the A are the squared masses of the exchanged quarks, the Ai and A2 are 



u 



the corresponding thresholds. It works efficiently in the case where the threshold 
parameters and s® differ from each other slightly. If we take the values s? = s° = 
(0.7 — 0.8) GeV 2 (in the case that the octet current J®(x) is chosen to interpolate 
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the 77 meson, see Ref. (25]) and s 2 = s® > M 2 o « 1 GeV 2 , the contributions from the 
a (980) are not taken into account properly, 



,0 „„0 

'i r 2 



i, / / Pqcd(si,s 2 )S(s 1 - gg) , , 

L«i - (p + 9) j [s 2 - <rj 

In the case of non-equal threshold parameters 7^ s°, we can take so = 
max(s ( j ) , s 2 ) with so small enough to avoid the contaminations from the high res- 
onances in either of the two channels, or take so = min(s ( j ) , s 2 ) with sq large enough 
to take into account the contributions from the ground states in either of the two 
channels. We have two choices in general, which can result in some uncertainties. 
In this article, we choose the current J^{x) to interpolate both the 77 and 7/ mesons 
to overcome the shortcoming, and take into account the contributions from the 7/ 
meson at the phenomenological side. 



3 Numerical result and discussion 

The input parameters of the light-cone distribution amplitudes are taken as A 3 = 0.0, 
f 3n = (0.45 ± 0.15) x 10~ 2 GeV 2 , u 3 = -1.5 ± 0.7, u; 4 = 0.2 ± 0.1, a x = 0.0, 
a 2 = 0.28 ± 0.08, a 4 = 0.0, r/ 4 = 10.0 ± 3.0 [SQ E3 M, 132], m u = m d = m q = (5.6 ± 
1.6) MeV, f n = 0.130 GeV, m w = 0.135 GeV, M n = 0.547 GeV, M v > = 0.958 GeV, 
M ao = 0.985 GeV, f v = 1.3f w , fa = 1.2/ w [40], and f ao = (0.21 ± 0.01) GeV [23]. 

The axial- vector current J M (x) has also non- vanishing couplings with both the 
pseudoscalar mesons 7/(1295), 7/(1405), 7/(1475), etc and the axial- vector mesons 
/i(1285), etc. The scalar current J(x) has also non- vanishing couplings with the 
scalar mesons cto(1450), etc. The masses and widths of those mesons are M v ^ 2 95) = 
(1294 ± 4) MeV, T v{1295) = (55 ± 5) MeV; M v{1405) = (1409.8 ± 2.5) MeV, r„ (14 05) = 
(51.1 ±3.4) MeV; M^ (1475) = (1476 ± 4) MeV, r^ (1475) = (87 ± 9) MeV; M /l(1285) = 
(1281.8 ± 0.6) MeV, r /l(1285) = (24.2 ± 1.1) MeV; M ao(1450) = (1474 ± 19) MeV and 
r ao (i45o) = (265 ± 13) MeV from the Particle Data Group [4"T] . 

From the experimental data, we can see that the ao channel permits larger thresh- 
old parameter than that of the 7/ channel. If we take the value so = max(s°, s° o ) = 
Sa ^ 1.7 GeV 2 , the contaminations from the 7/(1295) and /i(1285) are included in. 
We have to take the other choice, s = min(s°, s° Q ) = s° < 1.6 GeV 2 . It happens to 
be the ideal choice and reproduces the mass of the ao(980) with the conventional 
two-point QCD sum rules for the Borel parameter M 2 = (1.0 — 1.6) GeV 2 . 

In this article, we take the threshold parameter and Borel parameter as sq = 
(1.4 — 1.6) GeV 2 and M 2 = (1.0 — 1.6) GeV 2 to avoid the contaminations from the 
high resonances and continuum states as exp (—j^) = 0.2 — 0.4. In this region, 
the value of the coupling constant g is rather stable with variation of the Borel 
parameter, see Figs. (1-2). 
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Figure 1: The coupling constant g with the parameters M 2 and m u . 

In this article, we take the values of the coefficients a; of the twist-2 light-cone 
distribution amplitude <p n (u) from the conventional QCD sum rules [36], [39]. The 
(Ptt(u) has been analyzed with the light-cone QCD sum rules and (non-local con- 
densates) QCD sum rules confronting with the high precision CLEO data on the 
77* — ► 7T° transition form-factor [121 021 HU HSJ HH H7J . We also study the coupling 
constants ga oV n and g aori 'n° with the values a 2 = 0.29 and a 4 = —0.21 at /j, — 1 GeV, 
which are obtained via one-loop renormalization group equation for the central val- 
ues a 2 = 0.268 and a 4 = —0.186 at /i 2 = 1.35 GeV 2 from the (non-local condensates) 
QCD sum rules with improved model [17] . 

In the limit of large Borel parameter M 2 , the coupling constant g takes up the 
following behavior, 



It is not unexpected, the contributions from the two-particle twist-3 light-cone 
distribution amplitude ip p (u) are greatly enhanced by the large Borel parameter 
M 2 , (large) uncertainties of the relevant parameters presented in above equations 
have significant impact on the numerical results. The contribution from the two- 
particle twist-3 ipa( u o) is zero due to symmetry property. If we take the value 
m u = rrid = m q = (5.6 ± 1.6) MeV [39], the uncertainty comes from the m q is very 
large, about (33 — 64)%, and the predictive ability is poor, see Fig.l. From the Gell- 
Mann-Oakes-Renner relation, we can obtain f * m * = (0.027 ±0.003) GeV 3 [36], i.e. 

' m u +m d V ' ' 1 1 

m q w (5.6 ± 0.6) MeV, which can result in much smaller uncertainty. 

Taking into account all the uncertainties of the input parameters, finally we 




Qduo 



(15) 
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obtain the numerical values of the coupling constants, which are shown in Fig.2, 



9aor/'n 



for m q = (5.6 ± 1.6) MeV and 
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9aQr]ir° 
9aori'-K 



3.8tH GeV > 
3.llf?GeV, 
GeV, 

3.8l^GeV, 
3.118;? GeV, 
A.4+1 i GeV , 



(16) 



(17) 



for m q = (5.6 ± 0.6) MeV. The parameters of the twist-2 light-cone distribution 
amplitude <p n (u) obtained in Ref.[17j can change the value of the coupling constant 
slightly, less than 0.1%. 

In table 1, we list the values (not all) of the coupling constant g aori -K Q from different 
quark models and the experimental data. From the table, we can see that the values 
of the early estimations with the qq model, tetraquark model and KK molecule 
model deviate greatly from the experimental data [53j [Ml [55], we can't use them to 
identify the structures of the a (980) with confidence. Comparing with the values 
extracted from the radiative decay 0(1020) a (980)7 777^7 [5T| 1521 1531 1541 155] . 
the central value of our numerical result is reasonable and support the qq model. 



quark models and experimental data 


SWo(GeV) 


qq model [49J 


2.03 


tetraquark model [49J 


4.57 


KK molecule model [50], [8] 


1.74 


SND Collaboration [EH [52] 


3.11 


KLOE Collaboration [53, 54J 


3.0 ±0.2 


KLOE Collaboration [55J 


2.8 ±0.1 


light-cone sum rules (qq model) |25J 


2.6-3.4 


This work (qq model) 


1 +0.9 
°- i -0.7 



Table 1: The coupling constant g aoV n° from different quark models and experimental 
data. 



From the coupling constant g aorr K , we can obtain the decay width T 



Fao— >rpr° 



V 



P9aor)ir° 



127trf MeV for g = 3.8l?j GeV , 
127±|| MeV for g = 3.8±J;| GeV , 

Ml - (M v + m„Y] [Ml - (M v - m, 
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Figure 2: The coupling constant g with the parameter M 2 , A for m q = (5.6 ± 
1.6)MeV and B for m q = (5.6 ± 0.6) MeV. 

Comparing with the experimental data r ao (g 80 ) = (50 — 100) MeV from the Particle 
Data Group with very model dependent estimation (the decay width can be much 
larger) |41j, the central value and lower bound of our numerical result r^^^o = 
127^4g MeV are reasonable, however, the upper-bound is too large, we should reduce 
the uncertainties of the input parameters and m q (main uncertainties originate 
from them) before make definite conclusion. 

In this article, we take the point of view that the ao(980) is a scalar qq state. 
In Ref . |48j . the light nonet scalar mesons are taken as tetraquark states, and the 
coupling constants among the light scalar mesons and pseudoscalar mesons are cal- 
culated with the QCD sum rules. The numerical results indicate that the values of 
the coupling constants for the tetraquark states are always smaller than the corre- 
sponding ones for the qq states [22| [23]. 

The predictions listed in Table 1 are obtained from the phenomenological (po- 
tential) quark models [2111211501, and the resulting coupling constant g differs from 
the corresponding ones from the QCD sum rules greatly [2"2"j [231 HE]- Furthermore, 
those predictions also differ from the ones extracted from the experimental data 
significantly [5H [521 [531 EH [55] . In this article, we prefer the values from the QCD 
sum rules for consistence, i.e. if the nonet scalar mesons are tetraquark states, they 
have much smaller coupling constant g [221 E21 HE] • 

The scalar meson ao(980) may have small qq kernel of the typical qq meson 
size, strong coupling to the nearby KK threshold may result in some tetraquark 
components, whether the nucleon-like bound state or deuteron-like bound state. The 
tetraquark components may lead to smaller decay width, and smear the discrepancy 
between the (upper bound of) theoretical calculation and the experimental data. 
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4 Conclusion 



In this article, we take the point of view that the scalar meson ao(980) is a con- 
ventional qq state, and calculate the coupling constants g aoV n° an d g ao ri'iT with the 
light-cone QCD sum rules. Although the coupling constant g ao ri'iT has no direct 
phenomenological interest, we take into account the rj' meson to facilitate subtrac- 
tion of the continuum states to give more reliable sum rule. The central value of 
the coupling constant <? ao??7r o is consistent with the values extracted from the radia- 
tive decay 0(1020) — > ao (980)7 ~~ * rj-K '-/. The central value and lower bound of the 
decay width r^-^o = 127^4g MeV are compatible with the experimental data of 
the total decay width r ao ( 980 ) = (50 — 100) MeV from the Particle Data Group with 
very model dependent estimation (the decay width can be much larger), while the 
upper bound is too large. The scalar meson a (980) may have small qq kernel of 
the typical qq meson size, strong coupling to the nearby KK threshold may result 
in some tetraquark components, whether the nucleon-like bound state or deuteron- 
like bound state. The tetraquark components may lead to smaller decay width, 
and smear the discrepancy between the theoretical calculation and the experimental 
data. 
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Appendix 

We present some technical details in obtaining the spectral density at the phe- 
nomeno logical side, 

(v\J(o)\*°) = (^OlE/^^l^))^)!^ )!^)) 

= E/ (04 W)l g2 _ + ie H q ))(a( q )\ J(0)\AP)) 
= E/ ^(^KgK(p)) g2 _ ^ 2 + fe (a(g)l J(0)10) 

= E^/ (04^)1* / A^)l^») g2 _^ 2 + 2e 
= E/ (04 W) I / d 4 ygav(y)a(y)Ay)W<i)* (p)) M2 ie 

= E/(0i( 2 -w-^-^ ; /aMa 
- E 



(2tt)4 W ^ " ^M 2 -g 2 -* e 

gaf a M a 



M a 2 - (p' - p) 2 - it ' 

where we have used the completeness relation, 

d 3 q 



E/ (2^E |a(g))(a(g)l = 1 ' 



which corresponds to the normalization condition (a(q)\a(q')) = (27r) 3 2E5 3 (q — g'), 
the a's are the intermediate hadronic states with the same quantum numbers as the 
current operator J(0), the g a denote the corresponding coupling constant among 
the 7], a and 7r°, and (0| J(0)\a(q)) = f a M a . In the light-cone QCD sum rules, we 
often take the economical routine, 

(77|J(0)|7r°> = J2(v(p'M<i)Ap)) 2 _ l M2 . e (o(g)l^(o)|o> 

a Ha 

E gafaMg 
M 2 n -q 2 -ie' 

a a J. 

with a suitable definition (i](p')\a(q)ii (p)) = ig a - 
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The light-cone distribution amplitudes of the ir meson are defined as, 



{0\u(0)wd(x)\n(p)) 



\u(0)-f^g s G a(3 ( y vx)d(x)\n(p)) 



mix 



+iUml^- I' due- tup - x B{u) , 



(0\u(0)iKd(x)\n(p)) = 

{0\u(0)a^ 5 d(x)\ir{p)) = 

[0\u(0)(T a p'y 5 g 8 G IJH/ (yx)d(x)\n(p)) = 

(0\u(0)-f^ 5 g s G a p(vx)d(x)\Tr(p)) = 



2p ■ x J 

UK 1 



m u + m d J 



due~ lup - x if p (u) , 



due 



—iup-x 



Q(m u + m d ) J 

P a xp ~ Ppx, 



¥><r(«) , 



1 £ 2, 



p ■ x 

+Uml(ppg aii -p a g l 3 li ) 
Va i A ± (a i )e- ip < ad+vas) 

PaXp - PpX c 



p ■ x 

J Pa i V||(a i )e-^ (Qd+1,Q9) 



(19) 



where G a/ 3 = \e a p jXV G lxv and Vai = da>ida>2daz8{l — ct\ — a<i — ol%). 
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The light-cone distribution amplitudes are parameterized as 



ip p (u) 



V?37r(ai) 

V\\(ai) 
V\ (aA 



A{u) 



B{u) 



6u(l - u) |l + a x c} (2m - 1) + a 2 c| (2m - 1) + a 4 c| (2u - 1)} , 
l + {30r ?3 -^p 2 }c|(2 M -l) 

+ I -3^3 - -/ ~ -P 2 a 2 \ C!(2u - 1) , 



6u(l - u) <^ 1 



1 7 2 3 2 

1 



C|(2«-l) 



360a u a d a 9 jl + A 3 (a u - a d ) + ^ 3 -(7a 9 - 3) j , 

120a u a d a g (v 00 + v 10 {3a g - 1)) , 
120Q; M Q; rf Q; 9 aio(Q;d - a u ) , 

-30a 2 {h 00 (l - a g ) + h 01 [a g (\ - a g ) - Qa u a d ] 



+h 



10 



a 



/(! - <Xg) ~ ~ K + a2 d) 



30a g (a u - a d ) |/i o + h 01 a g + ^h w (5a g - 3)| , 

, f 16 24 20 
6u(l —«)< — + — a 2 + 20?73 + yr? 4 



+ 



+ 



117 10 

"15 + 16 " 27^ - 27* 



C|(2«-l) 



11 



210 



a 2 



135 



%^3 



C|(2m-1) 



18 
~5 



■a 2 + 217]^^ 



{2m 3 (10 - 15m + 6m 2 ) logM + 2m 3 (10 - 15m + 6m 2 ) log 
+uu{2 + 13mm)} , 

l + g 2 cl(2u-l)+g 4 cj(2u-l), 



u 



(20) 
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where 



3 

21 9 

20' 

21 

U 10 = -^^4^4 , 



^00 — w oo — ., • 



h 7 3 

7 3 
hio = -r}4V4 + 20 a2 ' 

18 20 
#2 = 1 + y a 2 + 60?73 + —r] A 



9 , , 

94 = -^a 2 -6r] 3 uj 3 , (21) 

113 3 

OHO* Q(O) C*2 2 (0 ancl Q(0 are Gegenbauer polynomials, r] 3 = lfL m ^+™d anc j 
p 2 = ( m -+ m ^ 2 [271 [2HI [291 EQl [3U EH [361 EZl EHl [39] . 

7T 

References 

[1] S. Godfray and J. Napolitano, Rev. Mod. Phys. 71 (1999) 1411. 

[2] F. E. Close and N. A. Tornqvist, J. Phys. G28 (2002) R249. 

[3] R. L. Jaffe, Phys. Rept. 409 (2005) 1. 

[4] C. Amsler and N. A. Tornqvist, Phys. Rept. 389 (2004) 61. 

[5] R. L. Jaffe, Phys. Rev. D15 (1977) 267, 281. 

[6] R. L. Jaffe, Phys. Rev. D17 (1978) 1444. 

[7] J. D. Weinstein and N. Isgur, Phys. Rev. D41 (1990) 2236. 

[8] J. D. Weinstein and N. Isgur, Phys. Rev. D27 (1983) 588. 

[9] N. A. Tornqvist, Z. Phys. C68 (1995) 647. 

[10] M. Boglione and M. R. Pennington, Phys. Rev. Lett. 79 (1997) 1998. 

[11] M. Boglione and M. R. Pennington, Phys. Rev. D65 (2002) 114010. 

[12] A. Bramon, G. Colangelo, P. J. Franzini and M. Greco, Phys. Lett. B287 (1992) 
263. 



15 



[13] P. J. Franzini, W. Kim and J. Lee-Franzini, Phys. Lett. B287 (1992) 259. 

[14] G. Colangelo and P. J. Franzini, Phys. Lett. B289 (1992) 189. 

[15] N. N. Achasov, V. V. Gubin and E. P. Solodov, Phys. Rev. D55 (1997) 2672. 

[16] A. Bramon, A. Grau and G. Pancheri, Phys. Lett. B289 (1992) 97. 

[17] A. Bramon, R. Escribano, J. L. Lucio M, M. Napsuciale and G. Pancheri, Eur. 
Phys. J. C26 (2002) 253. 

[18] E. Marco, S. Hirenzaki, E. Oset and H. Toki, Phys. Lett. B470 (1999) 20. 

[19] Z. G. Wang and S. L. Wan, Phys. Rev. D73 (2006) 094020. 

[20] Z. G. Wang, J. Phys. G34 (2007) 753. 

[21] Z. G. Wang and S. L. Wan, Phys. Rev. D74 (2006) 014017. 

[22] P. Colangelo and F. D. Fazio, Phys. Lett. B559 (2003) 49. 

[23] Z. G. Wang, W. M. Yang and S. L. Wan, Eur. Phys. J. C37 (2004) 223. 

[24] Z. G. Wang, Phys. Rev. D77 (2008) 054024. 

[25] A. Gokalp, Y. Sarac and O. Yilmaz, Mod. Phys. Lett. A19 (2004) 3011. 

[26] A. Gokalp, Y. Sarac and O. Yilmaz, Phys. Lett. B609 (2005) 291. 

[27] I. I. Balitsky, V. M. Braun and A. V. Kolesnichenko, Nucl. Phys. B312 (1989) 
509. 

[28] V. L. Chernyak and I. R. Zhitnitsky, Nucl. Phys. B345 (1990) 137. 

[29] V. L. Chernyak and A. R. Zhitnitsky, Phys. Rept. 112 (1984) 173. 

[30] V. M. Braun and I. E. Filyanov, Z. Phys. C44 (1989) 157. 

[31] V. M. Braun and I. E. Filyanov, Z. Phys. C48 (1990) 239. 

[32] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B147 (1979) 
385. 

[33] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B147 (1979) 
448. 

[34] L. J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rept. 127 (1985) I. 

[35] V. M. Belyaev, V. M. Braun, A. Khodjamirian and R. Riickl, Phys. Rev. D51 
(1995) 6177. 

16 



[36] P. Ball, JHEP 9901 (1999) 010. 

[37] P. Ball and R. Zwicky, Phys. Lett. B633 (2006) 289. 

[38] P. Ball and R. Zwicky, JHEP 0602 (2006) 034. 

[39] P. Ball, V. M. Braun and A. Lenz, JHEP 0605 (2006) 004. 

[40] L. Burakovsky and J. T. Goldman, Phys. Lett. B427 (1998) 361. 

[41] W. M. Yao et al, J. Phys. G33 (2006) 1. 

[42] A. Schmedding and O. I. Yakovlev, Phys. Rev. D62 (2000) 116002. 

[43] A. P. Bakulev, S. V. Mikhailov and N. G. Stefanis, Phys. Lett. B508 (2001) 
279. 

[44] A. P. Bakulev, S. V. Mikhailov and N. G. Stefanis, Phys. Rev. D67 (2003) 
074012. 

[45] A. P. Bakulev, S. V. Mikhailov and N. G. Stefanis, Phys. Lett. B578 (2004) 
91. 

[46] A. P. Bakulev and A. V. Pimikov, Acta. Phys. Polon. B37 (2006) 3627. 



[47] A. P. Bakulev, hep-ph/0611139 



[48] T. V. Brito, F. S. Navarra, M. Nielsen and M. E. Bracco, Phys. Lett. B608 
(2005) 69. 

[49] N. N. Achasov and V. N. Ivanchenko, Nucl. Phys. B315 (1989) 465. 

[50] N. N. Achasov and V. V. Gubin, Phys. Rev. D56 (1997) 4084. 

[51] M. N. Achasov et al, Phys. Lett. B485 (2000) 349. 

[52] M. N. Achasov et al, Phys. Lett. B479 (2000) 53. 

[53] A. Aloisio et al, Phys. Lett. B536 (2002) 209. 

[54] A. Aloisio et al, Phys. Lett. B537 (2002) 21. 

[55] F. Ambrosino et al, larXiv:0707.4609l 



17 



